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Abstract 

We investigate self-gravitating rotating solutions in the Einstein-Skyrme theory. 
These solutions are globally regular and asymptotically flat. We present a new kind 
of solutions with zero baryon number, which possess neither a flat limit nor a static 
limit. 



1 Introduction 

In non-Abelian field theories coupled to gravity particle-like solutions as well as black holes 
arise [1]. The latter are of importance as counterexamples to the no-hair conjecture. In 
recent years, in particular, various stationary rotating non-Abelian black holes have been 
studied [2] . However, the construction of stationary rotating particle-like solutions is still 
a difficult task. Although existence of such solutions in non-Abelian gauge field theories 
has been restricted [3, 4, 5], they can exist in the topologically trivial sector [6, 7]. On the 
other hand, stationary rotating soliton solutions in flat space have been obtained in the 
Skyrme model [8] and the U(l) gauged Skyrme model [9] in the nontrivial sector. 

The Skyrme model is a nonlinear chiral field theory in which baryons and nuclei are 
described in terms of solitons (so-called Skyrmions). Due to the long-standing difficulties in 
finding a satisfactory theoretical model for the interaction of baryons, much effort has been 
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devoted to the study of classical and quantised interactions of Skyrmions. However, the 
quantization of the Skyrme model is not only difficult since it is a non-renormalizable field 
theory; but also spinning Skyrmions must be considered which means that the solutions 
must consist of massive pions. This follows from the fact that the Skyrmion can only spin 
at a frequency up-to the pion mass before it begins to radiate pions [10]. Recently in [8], 
it was shown numerically that a good description of protons and neutrons can be achieved 
with spinning Skyrmions, provided the pion mass is chosen twice the experimental value. 

The rotating Skyrmion solutions of [8] are expected to persist, when the coupling to 
gravity is turned on gradually, analogous to the static Skyrmion solutions [11]. In the static 
limit, a branch of gravitating Skyrmions emerges from the flat space Skyrmion, when the 
coupling to gravity is increased from zero [11]. This branch terminates at a maximal value 
of the coupling parameter, when the coupling to gravity becomes too large for solutions to 
persist. A second branch of solutions exists which merges with the first one at the maximal 
value of the coupling parameter and extends back to zero. The solutions on the second 
branch possess a larger mass and they are unstable [12]. In the limit of vanishing coupling 
the solutions shrink to zero size and their mass diverges. As shown in [13], in this limit 
the Skyrmion solutions approach the lowest mass Bartnik-McKinnon (BM) solution of the 
SU(2) Einstein- Yang-Mills theory [14]. 

In this paper we investigate the stationary rotating generalisation of the static grav- 
itating Skyrmions. We show that in contrast to the static case, additional branches of 
solutions arise, which are not related to the flat space Skyrmion or to the BM solution. 
Most interestingly, we find a new kind of solution with zero baryon number, which exists for 
arbitrary (finite) coupling parameter, but does not possess a flat space limit. In particular, 
Section 2 presents the Einstein-Skyrme Lagrangian and the ansatz for the Skryme field 
and the metric, which lead to stationary rotating Skyrmions. In Section 3 the numerical 
solutions are discussed, while the conclusions are given in Section 4. 



2 Einstein-Skyrme Theory 

The SU (2) Einstein-Skyrme Lagrangian reads 

£ = — %z + — Tr (Kn Kn + -*-Tr ([K., K v \ [K^ K u \) + ^Tr ( l L±IL _ A M) 
16ttG 4 v M ; 32e 2 2 \ 2 J 

and its action is given by 

S = J C^d 4 x. (2) 

Here R is the curvature scalar, G is the Newton constant, k and e are the Skyrme model 
coupling constants, m n is the pion mass, and g corresponds to the determinant of the 
metric. The 577(2) Skyrme field U enters via = d^UU -1 . 

Variation of (2) with respect to the metric g^ v leads to the Einstein equations 

1 

G^ u = R^ v — —g^R (3) 
= 8-n-GT^, 
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where the stress-energy tensor is given by 

Tfiu = -yTr(V^ - l - 9lxv K a K^ ~^{g aP [K„ K a ] [K v , K^-g^ [K a , K p \ [K a , 

For stationary rotating solutions, two commuting Killing vector fields are imposed on 
the space-time: £ = d t and i] = d v in a system of adapted coordinates (t, r, 9, ip). In these 
coordinates the metric can be expressed in Lewis-Papapetrou form 

ds 2 = -fdt 2 + j (dr 2 + r 2 d0 2 ) + lr 2 sin 2 9 (dip - ^d?j * , (5) 

where f,m,l and u> are functions of r and 6 only. 

Then, the total mass and angular momentum are defined by 

respectively. Here S denotes an asymptotically flat hyper-surface, dV is the natural volume 
element on S, F is normal to S and = — 1. 

In order for finite energy configurations to exist the Skyrme field must tend to a constant 
matrix at spatial infinity: U — > 1 as — > oo. This effectively compactifies the three- 
dimensional Euclidean space into S 3 and implies that the Skyrme fields can be considered 
as maps from S 3 into SU(2). As the third homotopy class of SU(N) is Z, every field 
configuration is characterized by a topologically invariant integer B, which can be obtained 
as 

B = jB»k ll dV, (7) 

where is the topological current 

B» = J_ 2 ^ 2 e^Tr (K u K a Kp) . (8) 

This winding number classifies the solitonic sectors in the model and may be identified 
with the baryon number of the field configuration. 
For spinning Skyrmion the ansatz is of the form 1 : 

U — ni 1 + in 3 r z + in 2 (r x cos(</? + uj s t) + r y sm((p + u s t)) , (9) 

where r x ,r y , t z are the Pauli matrices; rii are functions of r and 6 only, satisfying the 
constraint C :— (1 — = 0, and the constant u s corresponds to the spinning frequency 

of the Skyrmions. 

1 Strictly speaking, the ansatz is neither stationary nor axially symmteric, since it depends explicitly on 
time and the azimuthal angle. However, the stress-energy tensor does possess the corresponding symme- 
tries. See also Ref. [4]. 
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When deriving the partial differential equations (PDEs) for the Skyrmion functions 
Hi, we have to take into account the constraint C — 0. This can be achieved by adding 
the constraint multiplied by some constant, say Co, to the Lagrangian and deriving the 
variational equations: 



„ dCJ^g dCJ^g dCJ^g 

dr^777{ T + Oe^77{ v 7{ 1" ^Tli 



d(d r rii) d{d e rii) 



dru 



-9 = 0. 



Then, the constant c can be obtained from the linear superposition Y, n iEi 
substituted back in the PDEs of (10). 



(10) 
0, and 



3 Numerical Solutions 



3.1 Parameters and Boundary Conditions 

Introducing the dimensionless radial coordinate x = ner, the gravitational coupling param- 
eter a 2 = 4tiGk 2 , the spinning frequency Cj s = cu s /Ke, and the pion mass rh n = m n / yj n/e, 
action (2) becomes 



K 



1 



^ + ^Tr (K, K») + ^Tr ([K„ K v ] [K», K»\) + ^Tr f^±^ - l' 



-g d 4 x , 
11) 



while the Einstein equations read: G^ v = 2a 2 T flu . We also introduce the dimensionless 
mass M = A^/[47r/c/e] and angular momentum J = J /[4ii/e 2 }. That way, the solutions 
depend only on the parameters a, Cj s and m n . For convenience we will rename Cj s — > uj s . 
At the origin, the boundary conditions are 

ni(0) = -l, n 2 (0) = n 3 (0) = 0, d x f\ = 0, d x l\ = 0, d x m\ = 0, u(0) = 0, (12) 

while for large x, since the asymptotic value of the Skyrme field is the unit matrix and of 
the metric is the Minkowski metric, we get 



ni(oo) — > 1, 7ii(oo)|i =2 ,3 — > 0, /(oo) — > 1, Z(oo) — > 1, ra(oo) — > 1, cj(oo) 

On the z-axis (# = 0) the boundary conditions follow from regularity 

d e n 1 \ e=0 = 0, n 2 (9 = 0) = 0, <9 e n 3 | e=0 = 0, 

def\o=o = 0, <9 e Z| e=0 = 0, <9 e m| e=0 = 0, <9 e w| e=0 = 0, 

while in the xy-plane (8 = 7r/2) from reflection symmetry 

<9 e ni| e=7r/2 = 0, d e n 2 \e =7T /2 = 0, n 3 (6> = tt/2) = 0, 

def\e=n/2 = 0, <%/| e=7r / 2 = 0, d e m\e =n / 2 = 0, <9 e ^| e=7r /2 = 0. 



0. (13) 



(14) 



(15) 
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In what follows we will encounter two special cases, the Bartnik-McKinnon solution 
and the non-trivial solutions in the vacuum sector. 

The first, is obtained after re-scaling x = ax and taking the limit of vanishing a. In 
this limit, the solutions are equivalent to the Bartnik-McKinnon one with lowest mass: 



where the gauge potential of the SU(2) Einstein- Yang-Mills theory is parametrized from 
w(x) via the relation A? — (1 — w(x))e ia jXj/(2x 2 ). 
The second is obtained by setting 



where the function h depends on x and 6. Regularity and finite energy of the solutions 
require that h vanishes on the z-axis and at infinity. We will refer to these solutions as 
'pion cloud'. 

3.2 Numerical Results 

The solutions are constructed using the software package CADSOL [15] based on the 
Newton-Raphson algorithm. In order to map the infinite range of the radial variable x 
to the finite interval [0,1] we introduce the compactified radial variable x — x/(l + x). 
Typical grids contain 70 x 50 points. The estimated relative errors are approximately 
~ 0.1%, except close to a max where they become as large as 1%. 

In particular, gravitating Skyrmions are constructed and their dependence on the cou- 
pling parameter a and the spinning frequency u s are studied for fixed pion mass: = 1. 
This sets a limit to the range of the spinning frequency uj s < 1. We start the discussion by 
a qualitative description of the dependence of the solutions on the gravitational parameter 
a for fixed spinning frequency uo s . Different branches of solutions exist which are char- 
acterized by their limit as a tends to zero. First, there are branches of solutions which 
tend to the flat space Skyrmions and to the scaled BM solution which we call (for obvious 
reasons) 'Skyrmion' and 'BM' branches, respectively. 

Second, branches of solutions exist which form a 'pion cloud' for large x as a — > 0. In 
this limit the solutions resemble a superposition of the 'pion cloud' with either a Skyrmion 
or the scaled BM solution and are called 'cloudy Skyrmion' and 'cloudy BM' branches, 
respectively. However, that for the 'cloudy' solutions the limit a — > is singular, as 
explained below. 

For fixed u s the branches exist up to a maximal value of a, where they merge with a 
branch of different type. In particular, the way different branches merge depends on the 
value of uj s relative to the critical value ujf ps 0.9607. So, for u s < ujf the 'Skyrmion' 
branches merge with the 'BM' ones; however for u s > uf, the 'Skyrmion' and the 'BM' 
branches merge with the 'cloudy Skyrmion' and the 'cloudy BM' ones, respectively. More- 
over, the 'cloudy Skyrmion' branches merge with the 'cloudy BM' branches only when 




(16) 



rii = cos(/i), n 2 = sin(/i), n 3 = 0, 



(17) 
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Next a quantitative description in terms of the dimensionless mass M and the value of 
the function l = 1(0) is presented. Since with vanishing a the mass diverges on the 'BM', 
the 'cloudy Skyrmions' and the 'cloudy BM' branches we also consider the scaled masses 
Ma and Ma 3 . 

Fig. 1(a) presents the mass M for the 'Skyrmion' branches (solid) which merge either 
with the 'BM' branches (dashed) or the 'cloudy Skyrmion' branches (dotted) when u s < uf 
and uj s > uj^ r , respectively. Note that, as a increases the mass decreases on the 'Skyrmion' 
branches, for small a; but diverges on the 'BM' and the 'cloudy Skyrmion' branches as 
a tends to zero. Also Fig. 1(a) shows the mass of the 'BM' branches merging with the 
'cloudy BM' branches (dash-dotted) when lo s > uof] and the mass of the 'cloudy Skyrmion' 
branches merging with the 'cloudy BM' branches when oj s < cuf '. 

The scaled mass Ma is plotted in Fig. 1(b). Note that, as a decreases along the 
'BM' branches the scaled mass Ma tends to a finite value which is equal to the mass 
of the Bartnik-McKinnon solution. Also, Fig. 1(b) shows the scaled mass aM of the 
'BM' branches as it merges with the 'cloudy BM' branches, and of the 'cloudy Skyrmion' 
branches as it merges with the 'cloudy BM' branches. Clearly aM diverges on the 'cloudy 
Skyrmion' and the 'cloudy BM' branches for vanishing a. 

Fig. 1(c) reveals that in the limit a — > the scaled mass Ma 3 of the 'cloudy Skyrmion' 
and the 'cloudy BM' branch tends to a unique value which depends only on uj s . 




Figure 1: The dimensionless mass M (a), the scaled masses Ma (b) and Ma 3 (c), and the 
value of / at the origin (d) as function of a 2 for several values of uj s . 
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Next we study the quantity lo of Fig. 1(d) in order to have a better understanding. 
Following a solution along a 'Skyrmion' branch which merges with a 'BM' branch, l 
decreases monotonically first along the 'Skyrmion' branch as a increases and then along the 
'BM' branch as a decreases, to take the value of the Bartnik-McKinnon solution as a — > 0. 
In contrast, when a 'Skyrmion' branch merges with a 'cloudy Skyrmion' branch, l reaches 
a minimum on the 'Skyrmion' branch and increases on the 'cloudy Skyrmion' branch as a 
decreases. On the other hand, when the 'BM' branch merges with a 'cloudy BM' branch, / 
increases with increasing a along the 'BM' branch until it reaches a maximum and decreases 
with decreasing a along the 'cloudy BM' branch. Finally, when a 'cloudy Skyrmion' branch 
merges with a 'cloudy BM' branch, l decreases monotonically if one follows the solutions 
first on the 'cloudy Skyrmion' branch with decreasing a and then on the 'cloudy BM' 
branch with increasing a. 

Also, while we observe that the scaled mass Ma 3 of the 'cloudy Skyrmion' and the 
'cloudy BM' branches tends to the same value as a tends to zero - this is not true for the 
quantity l . Thus, we conclude that the solutions approach different limits, though with 
the same (scaled) mass. 
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Figure 2: The dimensionless mass M and angular momentum J (a), and the scaled mass 
Ma 3 and angular momentum Ja 3 (b) as function of a 2 for lo s = 0.95 and uo s = 0.9608. 

In the following we compare the mass and the angular momentum. We here restrict to 
u s = 0.95 < oof and u s = 0.9608 > uf, as examples. In particular, Fig. 2(a) shows that 
when oj s = 0.95 < the angular momentum decreases monotonically on the 'Skyrmion' 
branch as a increases and on the 'BM' branch as a decreases, while it tends to zero on the 
'BM' branch as a — > 0. In contrast, for uo s = 0.9608 > a;" the angular momentum increases 
monotonically on the 'Skyrmion' branch as a increases and on the 'cloudy Skyrmion' branch 
as a decreases, while in the limit a — > the angular momentum diverges like a 3 . Fig. 2(b) 
presents the scaled mass Ma 3 and angular momentum Ja 3 for the 'BM' and 'cloudy BM' 
branches when u s = 0.9608 and for the 'cloudy Skyrmion' and 'cloudy BM' branches when 
uj s = 0.95. As a — > the (scaled) angular momentum tends to zero on the 'BM' branch, 
but takes finite values on the 'cloudy' branches. 

The new interesting feature of the rotating gravitating Skyrmions is the formation of 
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Figure 3: The functions n 2 (left) and n 3 (right) are plotted for 9 = 0, 7r/4, tt/2 for the four 
solutions with parameter values uj s = 0.95 and a = 0.1. 
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Figure 4: Analogous as Figure 3 for uo s = 0.9608. 
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the 'pion cloud'. This can be demonstrated by plotting the Skyrmion functions n 2 and n 3 
when a — 0.1 for uo s = 0.95 < uof and oo s = 0.9608 > uof as presented in Fig. 3 and Fig. 4 
(respectively). 

In particular, Figs. 3(a) and (b) show that on the 'Skyrmion' branch the solution is 
close to the flat space Skyrmion whereas on the 'BM' branch it is close to the (scaled) 
Bartnik-McKinnon solution. In contrast, the functions n 2 of the solutions on the 'cloudy 
Skyrmion' and 'cloudy BM' branch (presented in Figs. 3(c)) almost coincide at large x 
where the 'pion cloud' forms. However, for small x the functions n 2 are similar in shape 
to the corresponding ones of Fig. 3(a). Finally, by comparing the functions of the 
'Skyrmion' and 'BM' solutions (of Fig. 3(b)) with their 'cloudy' counterparts (of Fig. 3(d)) 
we observe that they are similar in shape for all x - which means that the 'pion cloud' 
does not reflect itself in the function n 2 . 

Fig. 4 presents the solutions on connected branches when oo s > cuf. Note that, the 
function n 2 of both the 'Skyrmion' and the 'cloudy Skyrmion' branch (plotted in Fig. 4(a)) 
almost coincide for small values of x and differ for larger x, where the 'pion cloud' is 
apparent. In contrast, the functions n 3 are almost identical for both solutions as shown in 
Fig. 4(b). Similar observations hold for the solutions of the 'BM' branch and the 'cloudy 
BM' branch plotted in Fig. 4(c)-(d). 

On the 'cloudy BM' branch the scaled BM solution in the core separates from the 
surrounding 'pion cloud' as a decreases. Therefore, in what follows, we show that pure 
'pion cloud' solutions can be constructed numerically by extracting the data of the 'pion 
cloud'. The ansatz for the 'pion cloud' solution is given by 

U = cos(/i) 1 + i sin(/i) (cos(y? + uj s t)r x + sin((/? + u s t)r y ) , (18) 

where the boundary conditions for the profile function h(x, 9) follow from finite mass and 
regularity requirements and read as 

h(0,6) = 0, h(oo,6)=0, h(x,6 = 0) = 0, d e h(x, 9 = tt/2) = 0. (19) 

For these solutions = so the chiral matrix can be regarded as a map from S 3 — > 
S 2 and thus, the 'pion cloud' solutions have zero baryon number. In addition, due to 
the boundary conditions of the profile function h(x,9), the 'pion cloud' solutions can be 
deformed continuously to the vacuum. 

Fig. 5 presents the scaled mass Ma 3 of the ' pion cloud' solutions and of the Skyrmions 
on the 'cloudy' branches as function of a 2 for fixed oj s = 0.95. Note that, for small a the 
masses of the 'cloudy Skyrmion' and 'cloudy BM' solutions coincide with the mass of the 
'pion cloud' solution. However, whereas the branches of the Skyrmion solutions exist only 
up to a maximal value of a, the 'pion cloud' solutions exist for arbitrarily large a. 

Indeed, as a increases the magnitude of h decreases linearly like I /a. So by expanding 
the Lagrangian up to quadratic order in h yields 



r A _I 

h 2a 2 2 



d x h d x h + \d e h d e h + - J. 2 h 2 - \ L s - -Y h 2 + m\h 2 
x 2 lx l sin f \ x J 



(20) 
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Figure 5: The scaled mass Ma 3 for the 'pion cloud' solutions and the Skyrmion solutions 
on the 'cloudy Skyrmion' and the 'cloudy BM' branches for u s = 0.95. 



which is equivalent to the Lagrangian of the rotating boson star 



'bs = — ~ ~y v (d^*d v $ + d v &dfi) - V{\$\), (21) 



B_ 1 

for $ = h(x, 9) e^-^st) an( j \/(|$|) = |$| 2 m 2 /2. The limit a — > 00 of rotating boson stars 
has been studied in Ref. [16] (though with different notation) where it was shown that the 
field equations become independent of the coupling parameter a after re-scaling h — h/a. 
Moreover, it was argued in [16] that several branches of solutions exist in certain ranges 
of u s which suggests that several branches of 'pion cloud' solutions might (also) exist, at 
least for large values of a. 

In the limit a — > the scaled mass Ma 3 takes finite values. By introducing the scaled 
radial coordinate £ = ax, the Skyrme field equations reduces to the constraint 

sin 2 (/i) [cos(/iV 2 - fml] = , (22) 

as a — > 0. This implies that the Skyrme field function is either h = or cos(/i) = /m 2 /^ 2 . 
Note that none of these solutions can hold globally; the former one yields the trivial solution 
and the latter is not consistent with the asymptotic boundary conditions h — > and / — > 1 
for u 2 s < m 2 . However, both can hold locally. Indeed, we find the solution 

2 

Tfl 

cos(h) = {^}GD = (0,£o]x(0,7r) 

^ s 

h = , elsewhere , 

where £0 depends on u; s . Although this solution for h is not continuous at the origin and 
on the z-axis, the metric functions are continuous globally. Moreover, in the domain V the 
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function h depends only on the radial coordinate £. Consequently, the metric is spherically 
symmetric and, for £ > £ 0) given by the Schwarzschild solution. 

Thus, in the limit a — > the 'pion cloud' becomes confined to the finite domain T>. 
Outside of D the metric is the vacuum solution determined by the connecting conditions at 
£o- To see the physical picture, however, we have to return to unsealed coordinates. Then 
the domain V extends over the whole space, except the z-axis and infinity. Consequently, 
the 'pion cloud' occupies an increasing volume in space as a decreases. 

The profile function h is plotted in Fig. 6(a) as function of the cylindrical coordinates 
p = xsin9, z = xcos9 for oj s = 0.95 and a = 0.07 and in Fig. 6(b) as function of the 
scaled coordinates p' = £sin6>, z' — £ cos# for such parameter values that the constraint is 
almost satisfied (ie. oo s = 0.95 and a = 10~ 4 ). 



h(p,z) 



(a) 




(b) 




Figure 6: (a) The 'pion cloud' profile function h(p, z) when uj s = 0.95 and a = 0.07. (b) 
Same as (a) for a = 10~ 4 in scaled coordinates p' — £ sin 0, z' — £ cos 9. 



Finally we state that the 'pion cloud' solutions do not exist for arbitrarily small u s . In 
fact, we observed that the coefficient of the second order derivative term in the Skyrme 
field equation 



^) sm 2 (h)-f)r 2 lsm 2 9 



9 z rr h + ^d' z ee h )+••• = (23) 



develops a zero at some point on the 9 = it/ 2 axis, when uj s decreases to a critical value, 
and no solution exists for uj s below the critical value. Therefore, the 'pion cloud' solutions 
do not possess a static limit. 



4 Conclusions 

We studied stationary rotating solutions of the Einstein- Skyrme theory. These solutions 
are asymptotically flat, globally regular and axially symmetric. Branches of stationary 
rotating Skyrmions emerge from corresponding branches of static Skyrmions, when the 
rotational frequency is increased form zero. If the rotational parameter oj s is smaller than 
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a critical value, the rotating Skyrmions on these branches behave similar to their static 
counterparts, when the coupling to gravity is varied. However, additional branches of 
solutions exist, which do not have a flat space limit. These branches are characterised by 
the formation of a 'pion cloud' for small values of the coupling a. 

In addition to the rotating Skyrmions with baryon number one, we found new solutions 
in the topologically trivial sector. These 'pion cloud' solutions are also asymptotically flat 
and globally regular, but possess neither a flat limit nor a static limit. In contrast to the 
Skyrmions the 'pion cloud' solutions exist for arbitrary coupling, ie. < a < oo. 

Axially symmetric rotating Skyrmions with higher baryon number B > 1 should easily 
be obtained by the replacement (p — > B(p in the ansatz Eq. (9). 

Einstein-Skyrme theory also possesses black hole solutions. So far the Skyrmion black 
holes have been studied in the static limit only. Families of black hole solutions emerge 
from the globally static regular Skyrmions, when the horizon radius is increased from zero. 
Similarly, we expect several families of stationary rotating Skrymion black holes to emerge 
from the globally regular rotating Skyrmions on the different branches obtained here [17]. 
Moreover, one may speculate that also stationary rotating black holes emerge from the 
'pion cloud' solutions. 
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